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ABSTRACT 

There is a need for this country to be able to assure itself that 
an international ban on nuclear testing can be made effective. This 
thesis studies one small portion of that problem - to wit, the procedure 
for classifying an underwater detonation which has vented and which has 
been detected by a hydroacoustic network. 

A MEASURE OF EFFECTIVENESS is developed, and after a study of some 
parameters such as height of search, search spacing, search width and 
current movements optimizing search plans are deduced. 

One of the major difficulties in such a quantitative analysis is 
the non-deterministic nature of many of the parameters. The ideal 
search spacing was found to depend directly on the surface radioactive 
pool diameter. In order to develop plans this essential factor must be 
known. Unfortunately this depends on many variables over which the 
searcher has no control or knowledge. An averaging procedure is adopted 
and a factor established to compare spacings and determine which meets 
the criterion dictated by the MEASURE OF EFFECTIVENESS, : 

The writer wishes to express his appreciation for the assistance 
given him by Professor H. E. Handler, Professor R. Read and Professor 
E. M. Ward; for the assistance and encouragement given him by Professor 
W. P. Cunningham; and particularly to Professor E. J. Stewart for the 
kind and charitable way in which he made his time and advice available. 
Lastly, the writer wishes to thank Mrs. Martha Olson of U. 5. Naval 
Radiological Defense laboratory for making available hydrographic 


information compiled by her. 
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1. Introduction. 


In October 1960 a technical symposium was held by the Advanced 
Research Projects Agency (ARPA) acting in behalf of the Department of 
Defense for the purpose of informing those in attendance of the scope, 
status and planning underway for Project VELA. Project VELA consists of 
research and development that is aimed at the improvement of our capa- 
bility of detecting and identifying nuclear detonations that are con- 
tained below the surface of the earth, or underwater or that occur in 
space above the atmosphere. 

It is well recognized that it is possible for a clever and deter- 
mined nation to conduct subsurface and high altitude nuclear weapons 
tests without a significant risk of being detected. Detection of 
nuclear detonations on the surface or within the atmosphere is consid- 
ered well in hand and not included in VELA. Ever since the Conference 
of Experts held in 1958 in Geneva, international groups have been 
attempting to establish methods and networke which would assure the 
participating nations that agreements to ban nuclear testing could be 
effectively enforced and thereby be the first step in solving one of the 
most important problems of our time - ARMS CONTROL. 

The hope of Project VELA is to improve our detection and classifi- 
cation capability by at least one order of magnituwie, and perhaps by as 
much as three, over a three year span. Such a quantum jump is considered 
feasible if sufficient independent scientific interest can be aroused 
while assigning a high priority to the investigation in as many of the 
govermment laboratories 7 have some measure to contribute. 

The Project was subdivided into three major categories: 
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VELA HOTEL - Detection of very high altitude nuclear detonations by 
means of detectors mounted in satellites; 

VELA SIERRA - Development of techniques and equipment to detect and 
identify high altitude nuclear detonations by ground based instmmenta- 
tion; 

VELA UNIFORM = Improve, markedly, equipment and techniques for de- 
tecting and identifying subsurface nuclear explosions from a distance as 
well as to develop rapid, economical techniques for conducting on-site 
inspections to determine whether a suspicious event was actually a nuclear 
explosion or not. 

Project VELA UNIFORM includes a study of underwater detonations as 
part of its "subsurface" objective. The Naval Electronics Laboratory 
(NEL) at San Diego was assigned the task of investigating the detection 
of underwater explosions and the Naval Radiological Defense Laboratory 
(NRDL) was contracted to devise a searching procedure for obtaining san- 
ples to classify explosions detected by a network devised by NEL. 

In a CONFIDENTIAL report /6/ NEL outlined a distribution of hydro- 
acoustic stations which would be capable of detecting underwater blasts 
in most areas of the world. However, this network has practically no 
ability to classify explosions as nuclear or non-nuclear. 

NRDL has pursued the problem of classifying a reported detonation 
by locating the scene by means of an appropriate search based on the fix 
information received from hydroacoustic and/or other correlated equip- 
ment, obtaining samples if possible and analyzing pertinent character- 
istics of the area. 

The classification problem splits into two distinct categories re- 
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flecting two completely different physical characteristics of the scene. 
In the first case, the underwater nuclear explosion vents to the surface 
and scatters the attendant radioactive debris on the surface and probably 
into the air. In the second case the radiation is contained below the 
surface due to a combination of depth of explosion and yleld of the 
device. This thesis is devoted to analyzing the first of these two 
categories. 

2. Statement of the problem. 

The objective of this study is to determine the optimum searching 
procedure for locating the site of an underwater explosion which has 
vented, in order to classify it as nuclear or non-nuclear. 

3. Means of identification. 

The most positive means of confirming that an explosion is nuclear 
is by a radiochemical analysis of a sample of the affected area. 

After an explosion has been detected and a radiating sample obtained, 
it remains to distinguish the origin of the radiation. In spite of a 
radioactive sample a nuclear detonation could be disclaimed and the 
activity attributed to one of the following: 

a. Residue from previous explosions, 

b. Waste dumping, 

€. Reactor discharge, 

d. Accident. 

Radiochemistry can detect fission fragments and normally determine 
the time of detonation by extrapolating back from ratios of certain iso= 
topes present in the sample. This information correlated with seismic, 
hydroacoustic, air acoustic, and/or electromagnetic information of time 
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and location of detonation can readily refute a claim of residue from 
previous explosions or waste dumping as the source. 

Normal reactor discharge does not contain fission fragments and 
would easily be repudiated as a cause for contamination. However, leak- 
age from a reactor primary system could be a source of deceiving activity 
in the sample, but would not explain the detected explosion at the loca- 
tion. 

An explosion of a reactor in water would have attendant debris, but 
with the burden of proof to counter such a claim upon the investigator, 
failure to find debris is unlikely to be satisfactorily conclusive that 
such an explosion did not take place. 

Any surface sample is subject to prior masking by the dumping of 
fission fragments to confuse analysis. A determined effort to make 
radiochemical analysis of surface samples inconclusive is almost bound 
to be successful in a water contaminated area if the deceiver is willing 
to invest enough effort. Contaminated air samples from venting of under- 
water bubbles are less vulnerable. 

Underwater pulsing of the bubble on its route to the surface leaves 
a trail of discrete thin laminas of radioactive material. To reduce the 
possibility of confusion of the surface sample by a masking attempt a 
further effort should be made to sample some of these lower strata. 

sampling procedures should recognize biological phenomena which can 
corroborate other evidence. Certain algae have a known affinity for 
Iodine 131 which has a short half life and can be identified as a fission 
fragment. Plankton can also produce evidence of an event. If there are 
land and shallow water in the vicinity, molluscs and tropical plant sam- 
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ples Eould predues valuable information. 

Of course, specifically, well equipped search vehicles could arrive 
at reliable conclusions at the scene with pulse height analysis equip- 
ment, air drop water probes and well schooled personnel. Howewer, it is 
not realistic to imagine planes particularly outfitted for this searching 
chore. Our forces must be more flexible for the large variety of mis- 
sions they may be called on to perform. It is envisioned that patrol 
craft could be rigged with the detector described in USNRDI-TR-525/1/ or 
that bases from which these planes operate could have these meters avail- 
able. However the more sophisticated equipment is not likely to be as 
well distributed. 

Actual meter readings would be of some assistance in correlating 
with the time of the explosion by determining decay rates; however, con- 
tinued dilution due to wind and sea movements will make analysis of this 
data much less reliable than water samples. 

4. Detection capability. 

It is worthwhile to consider some of the major factors involved in 
the transmission of sound waves in water. 

The sound reaching a receiver experiences three kinds of degradation: 

a. Divergence - For a point source in a homogeneous medium the 
energy is dispersed in a sphere as it is radiated out so that a receiver 
perceives only an amount of energy proportional to the percent of area 
of the expanding sphere that it can observe. 

b. Attenuation - There is a damping of the energy because the pres- 
sure is not in phase with the condensation due to viscosity, heat conduc- 
tion and molecular absorption (due to magnesium sulphate). 

c¢. Anomalies - There are many anomalies which contribute to the 
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peculiarities of underwater sound transmission but the principal ones are 
reflection and refraction. 


The velocity of sound in a fluid (C) = 8} 
| 


where B = Bulk Modulus 


P = Density 


As Temperature goes up P decreases; 

As Salinity goes up B increases; 

As Pressure increases P increases but B increases more. 

Therefore the velocity increases as temperature, salinity and pressure 
increase. 

Pressure increases with depth but temperature noreally decreases, and 
temperature reductions have an overriding effect on the velocity. In the 
Atlantic a temperature of 2°C is reached at about 4000 feet and remains 
fairly constant from then on. Salinity remains relatively constant and 
velocity of sound increases with depth due to the pressure increase. In 
the Pacific the same conditions obtain at about 2000 feet. Near the 
equator this condition of constancy is reached at a lesser depth, 

Near the surface radical temperature changes are not at all infre- 
quent. Surface disturbances of storms mix the water for a considerable 
depth and sometimes send warm water down while the colder water is 
brought to the surface. This produces temperature inversions which have 
a drastic effect on sound transmission. This will be discussed in more 
detail below. 

As the velocity increases the ray (perpendicular to the wave front) 
curves downward. Conversely, if the velocity decreases the ray curves 


upward. 
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Consider the. following ray diagrams where the left depicts the ten- 
perature changes with depth and the right shows the ray path (Go is 


original direction of propagation). 
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Fig. 1. Effect of temperature changes with depth on sound ray path. 

Sound channes occur when there is a minimum velocity at some depth 
with higher velocities above and below. The region of minimum velocity 
is called the axis. Sound originated on the axis and propagated along 
the axis will not diverge spherically but only cylindrically which means 
that the losses are cut down. Additionally, there is no major anomaly 
to degrade the wave. 
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b. Pseudo Sound Channel 
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Fig. 2. Sound channels. 

Sound sources do not have to be in the sound channel in order for 
sound to be propagated along the channel. Rays entering sound channels 
at small enough angles (usually less than 15°) with the horizontal may 
be entrapped and then act as if they originated on the axis. 

These few facts may help the reader to appreciate the following 
comment paraphrased from the preliminary NEL report /7/: 

Historically it is apparent that even small explosions under idea 

conditions can be detected at ranges of many thousands of miles , 

whereas, the most powerful explosions under other conditions may 

1 MILES will be NAUTICAL MILES (6087 feet) throughout this thesis. 
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go undetected. For example; one pound of TNT exploded on the sound 

axis in the open ocean would be expected to produce a stronger sig- 

nal than a 100 Megaton explosion under other extreme conditions. 

Obviously, then, the hydroacoustic detection of a nuclear burst is 

a complex problem and a generalized statement of the capability of 

such a system is not possible. 

The interfaces between the air and the water or the land and the 
water provide a poor match for transmission of sound waves (high coeffi- 
cient of reflection). As a result surface bursts or underground bursts 
are not well detected through hydroacoustic means just as underwater 
bursts are not well detected by seismic devices. For this same reason 
stations are not able to monitor inland seas and locations in lagoons or 
with seamounts between the event and the receiver. 

The velocity of sound in water is about 4800 feet a second and it 
is possible that some stations will receive a signal as long as one hour 
after detonation (less than 3000 miles). It will take at least two more 
hours to correlate station information, evaluate, and call an alert. 

There are many sources of confusion in detecting. Hydroacoustic 
stations will be unable to distinguish between chemical explosions and 
atomic explosions. Earthquakes and underwater volcanic disruptions when 
detected will be a concern and frequently non-discernible from man made 
explosions. The splash of missiles in deep water, ship propulsion noises 
and enviromental sounds characteristic of the ocean can be distinguished. 

In the Joint Committee Hearings /8/ Dr. H. Brown of Lawrence Radia- 
tion Laboratory estimated that there were 5000 continental earthquakes 
per year that could be suspected of being nuclear explosions of between 
one and five kilotons and that only a few percent could be identified as 
earthquakes by present criteria. About 100 disturbances would be sus- 


pected of being nuclear explosions of five kilotons or more. Since most 
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of the continental seismic false alarms are caused by these earthquakes 
it seems unlikely that the coe of hydroacoustic false alarms due to 
natural causes would be this order of magnitude. Earthquake disturbances 
are not likely to be transmitted with too much energy through the inter- 
face with the water. Volcanic disruptions below the water may be mis- 
leading, but these are not near as numerous as continental earthquakes. 
The principal cause of false alarm in the hydroacoustic network is most 
likely to be chemical explosions transmitted by the Sound Channel mecha- 
nism. 

5. Alternative courses of action. 

With our objective clear in mind - to determine whether a reported 
explosion was nuclear or non-nuclear - our problem is reduced to search- 
ing the area in question to determine if there is evidence of nuclear 
activity. Only by not obtaining evidence of nuclear activity can we 
assume that the disturbance was non-nuclear (ignoring the possibility of 
non-venting nuclear explosions for this study), and even then the conclu- 
sion is not absolute. The certainty with which a conclusion can be 
reached depends on the thoroughness of the search. 

The courses of action that are open revolve about the approach to 
the development of a search plan. 

a. Should the plan be based upon a random search approach, consid- 
ering the target to be equally likely in some large area prescribed by 
hydroacoustic fix information? This would permit parallel sweeps through 
the area in an orderly fashion by several aircraft with no danger of in- 
terfering with each other. It would assure a thorough coverage of the 
area because navigation is simple and reliable. 
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b. Should the plan attempt to optimize search effort with respect 
to time? This would permit target distribution to be any shape that the 
physical situation dictates. The pattern would be formed to search the 
areas of greatest probability continuously with the distributim of the 
fix being the "a priori" function and an "a posteriori" function formed 
after each search. This approach would be best for the following con- 
siderations: 

(1) Search craft may have to be recalled prior to completion of 
the plan or there may be a cost or utility consideration. 

(2) International ground rules may limit time to search a given 
area or it may be desired to spend as little time as possible in a given 
area due to political or other considerations. 

(3) There is the ever present marginal case where just enough 
radioactivity is present at one time in sufficient intensity (because of 
yield, depth of detonation, and/or dispersion) to be detected, but which 
will not be there at some later time. 

A disadvantage to this approach to the search plan is that, unless 
target location distribution is in fact uniform, the navigation problems 
reduce the reliability of conforming to a given pattern, and holidays 
are left in the coverage. The probability of detection may in fact be 
less than a random pattern, although the same coverage effort is invested, 
because of those holidays. Secondly, there is a peculiarity in the 
radioactive pool contour under many conditions which reduces the proba- 
bility of detection at certain times. The character of many pools is as 
follows: 

Upon initial venting a high intensity area is established for a 
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certain distance around “eround zero." After some time passes many of 
the fission fragments sink, and fringe area activity is diluted beyond 
recognition. This decreases the contour of a given low isodose line. 
Then the particles that are suspended near the center spread out with 
water and wind movements and the low isodose contours spread out again, 
even beyond the original area. 

It is impossible to fill in details of "how much?" or "how long?" 
or "what intensity?" because all answers depemd on the physical condi- 
tions of venting and water movement. The point of mentioning it here is 
to acknowledge that it could affect the capability of our searching 
effort. 

6. Measure of effectiveness. é 
This situation is different from many military operations where the 
accomplishment of the objective is imperative, regardless of the cost. 
The search would be conducted during peacetime and there may be concern 
for the expense involved, parti¢ularly if many false alarms are Siiiuntenll 
There is always a possibility of search vehicles having to be recalled 
because of more urgent commitments, logistic considerations or forbidding 
weather. As previously mentioned, international agreements may limit 
searching time. 
Finally, there is a matter of using consistent logic in approaching 
the problem. As presented so far, and as the physical situatim dictates, 
the problem admits the probablistic nature of the fix information. Again 
ignoring the non-venting situation, there can never be an absolute assur~ 
ance that the explosion was non-nuclear, but the probability of this fact 
approaches unity as the search effort approaches an extremely large value. 
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Therefore, the search will terminate, if a pool is not located, at some 
time before the evidence is absolutely conclusive that no pool exists. 
This will always be the case and the amount of search effort, or synony- 
mously - time, is inherent in the problem and cannot be ignored. 

Under these considerations the MEASURE OF EFFECTIVENESS is to maxi- 
mize the probability of detection per time of search effort. 

7. Ildst of variables to be considered. 

a. Height of search. This is discussed in detail in Appendix I. 
The conclusion is to fly as low as operationally feasible. A plot of 
intensity versus height in Appendix I shows the loss of intensity for 
altitudes from 0 to 40,000 feet. 

be. Detector sensitivity. Appendix III is based on a detector which 
is capable of distinguishing 0.2 microroentgens at the meter. In order 
to meet this requirement the meter envisioned will be collimatedto reduce 
the effect of cosmic radiation. 

c. Search spacing. In Appendix II the search theory involved is 
discussed, and it is generally concluded that the spacing depends upon 
the search width (which is the effective visibility in the definite range 
law) and the standard deviation of the location errors. 

d. Search width. The radioactive pool on the surface of the water 
coupled with the width of the detector view determine the search width. 
This is covered in detail in Appendix III. 

e. Target motion. The movement of the radioactive pool, which is 
the object of the search, will depend principally on currents. Appendix 
IV investigates the variability of the water movement and concludes that 
the speed of movement is predictable but the direction is random. 
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g. Probability of detection. In Appendix III the probability of 
detection is computed for the ideal spacing - given the search width. 
Since the search width is not known in a realistic situation an aver- 
aging procedure must be resorted to. This is developed in paragraph 11 
below. 

h. Search aircraft characteristics. 

1. Navigation errors in locating the center of the fix. 

j. Time late at the scene. 

8. Hydroacoustic fix distribution. 

The location of an explosion will be estimated by correlation be- 
tween hydroacoustic detection stations. Their estimates will depend on 
physical condtions prevailing, but for this study the errors have been 
classed as having a standard deviation of 2.5 miles or 25 miles. The 
purpose of assuming deviations an order of magnitwie apart is to see if 
the search plan is sensitive to such extremes. The distribution is as- 
sumed normal and the probability is represented by: 


f(x,y) = f(r) = she KPa} 
20" (az 


On the basis of this assumption Appendix II develops the ideal 
search spacing as 0.75 (EC)2 where E is the Effective Visibility and ¢ 
is the standard deviation. 

9. Navigation error. 

Where there are no visual navigational aids available an aircraft 
can expect to arrive at a given location in accordance with some proba- 
bility distribution. According to Bowditch /4/ celestial errors in air 
navigation of five to ten miles are considered normal for favorable con- 


14 


ditions (pg 676 paragraph 2806), Electronic navigational aids, where 
available, yield more accurate positioning. In this model it is assumed 
that the error due to navigation is normally distributed with a standard 
deviation of five miles. The variance is additive to the hydroacoustic 
fix variance since they are both normal and independent. When the fix 
standard deviation is 2.5 miles the total location standard deviation is 
5.6 miles. When the fix deviation is 25 miles the location deviation is 
25.5 di lea. 

These figures demonstrate that the navigation error has a substan- 
tial effect on the problem when the fix error is small, but it is rather 
inconsequential when the fix error is large. 

10. Search patterns. 

Appendix IV discusses curreits in detail and for this study the set 
is considered equally likely in any direction at a given location and 
time. The drift, on the other hand, can be predicted reasonably well. 
Koopman /2/ in sections 1.6 and 7.3 discusses a non-uniform distribution 
of targets. When the fix deviation is 2.5 miles the situation falls 
closely into the category of a circular distribution of the target.with 
movement at a known speed in a random direction. The theory is to search 
on the peak of the distribution, then outside, then inside, ... etc. 

This pattern will give the maximm probability of detection per unit time 
for locating a moving target whose initial fix is not definite. See 
Appendix V for details of this plan. 

When the location deviation is of the order of 25 miles or more and 
the time late is less than two days the target can be treated as station- 
ary and the search pattern that applies is described by Koopmen /2/ in 
paragraph 7.3.2. A set of axpanding squares ~ squares of uniform cover- 
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age - are superimposed in the search pattern. See Appendix VI for de- 
tails of this plan. 
11. Seareh width. 

This one parameter remains to be determined. Appendix III discusses 
the geometry for determining the search width and it concludes that it is 
the width of a uniform pool of minimum intensity minus the detector view. 
Unfortunately, there is no way to determine the pool size. There are 
really not enough statistical data available to generally characterize 
pool contours for various depths, yields, sea and wind conditions an as 
a function of time. 

a. Large location error. 

As a matter of interest the following tables have been computed 
to consider the effect of using a search width for one size pool when 
the pool is actually some other size. These calculations are based on 


the estimated probability of detection given by Koopman /2/. 


Po 1 ~ {1 ta exp =nw 
s ) s 
n = Number of searches made (one for this table). 


w = Search width. 


Spacing of pattern. 


g 
As concluded in Appendix II the spacing is £75 (ECE and for this 


problem W = E. The standard deviation considered is 25.5 miles. 
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TABIE 1. Analysis of probability of detection for given search spacings. 


Search spacing (S) = 4.31 miles - Computed for a 2 mile diameter pool 


Pool Diameter Search Width W Probability of 


(Pool of uniform minimm intensity ) (W = E) S Detection 
2 Teo 302 0.036 
4 Se ~766 018 
6 D3 1.23 0348 
8 73 1.692 «505 
10 9.3 2616 0634 
20 19.3 4.48 939 
30 29.3 6.8 0991 
spacing (S) = 6,86 =~ Computed for a 4 mile diameter pool 
2 1.3 -1896 015 
4 aeo 481 085 
6 503 07725 aa 
8 763 1.062 0287 
10 9.3 1.352 394 
20 19.3 2.813 o7TT14 
30 29.3 4.27 ©9262 
Search spacing (S) = 8,7 - Computed for a 6 mile diameter pool 

2 “te 01495 01 

4 ae 03795 055 
6 aes «609 n> 
8 73 -839 2205 
10 9.3 1.068 2289 
20 19.3 222 0652 
30 Zao 30375 08513 
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Search spacing (S) = 10,2 - Computed for an 8 mile diameter pool 


Pool Diameter 


0 Aa hDNY YF NN 


1 
20 


30 


Search Width 


1.3 
303 
5 3 
73 
93 
19.3 
29 3 


W 
S 


lero 


324 
052 

2715 
aie 


lieeo> 
Zeon 


Probability of 
Detection 


0075 
04 
2095 
0161 
ee2a 
0564 
o78 


Search spacing (S) = 11.5 - Computed for a 10 mile diameter pool 


ao GF fF NWN 


10 
20 


30 


1.3 
33 
503 
73 
923 
19.3 
e780 


0113 
0287 
461 
0634 
~808 


1.678 
2455 


- 008 
034 
076 
0132 
0191 
9 

0123 


Search spacing: (S) = 16,6 - Computed for a 20 mile diameter pool 


oO Oma DBF YF N 


30 


1.3 
33 
303 
723 
93 
19.3 
29.3 


.0783 
01987 


0319 
Ak 
056 


1.163 
1.764 


«005 
00175 
004 

«973 
2107 
0325 
0525 





Search spacing (S) = 20,4 - Computed for a 30 mile diameter pool 


Pool Diameter Search Width W Probability of 
. S Detection 

2 1.3 0637 «002 
4 3 «3 01618 01 
6 53 026 029 
8 73 358 05 

10 9.3 456 078 

20 ie 946 e243 

30 29 3 1.436 oA21 


Figure 3 is a plot of the probability of detection for each of the 
search spacings under considerations against the widths of the pools. A 
plot against search widths would be the same shape displaced to the left 
O.7 miles. Figure 4 utilizes the same data but uses SPACING vice POOL 
DIAMETER as the abscissa. 

In order to consider this nebulous problem more quantitatively, 
assume that any search width from.0 to 30 miles is equally likely and 
investigate the mean probability of detection for each search spacing 


over the span of search widths. 


30 
= f p(w)f (widw = It 7. ) era 
QO 30 


P 

P= af freocy -W y expla 

_ 30 

P= 1 ist exp(-W) + S xp (i) war 1 
30 S) (S 1. 

pe 


a] 30 +5 exp (30) +S exp(-30) (30+1) ~ 23 | 
30 S) Ss) (8 


The Measure of Effectiveness is again called upon to decide which 
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spacing affords the best approach. It is desired to select the spacing 
which yields the maximum probability per unit time of search effort. The 
spacing (S) is inversely proportional to the time of search effort. 

Hence, by maximizing pS the criterion previously established will be 
satisfied. 

TABLE 2. Comparison of spacings using the MEASURE OF EFFECTIVENESS (LARGE 


LOCATION ERROR, 


Associated Spacing Ave. Prob. of Detect. Measure of 
Pool Diameter (SS) _ (5) Effectiveness (SS) 

2 4.31 (071355 3.07 

4 6.86 0561145 3.855 

6 8.7 "447058 4.10 

8 10.2 e409 a, fonlier 

10 11.5 23633 40175 

20 16.6 0239 32965 

30 20.4 183 3.73 


From the above table it can be seen that pS is maximized around a 
spacing of 11.5 miles. 
By differentiating the expression for pS with respect to S and 


solving for the maximum the following is obtained: 
d(ps)= 4 < S+S* oxp(-30)4$- exp(-30)(30+ )- | 


ds 30 30 )(s 15 
d (pS )=0=1+exp(-30)+_1 S exp(-30)+30 exp(-30)+2 exp(-30)+_S exp(-30) 
as . S ) 15 Ss )s S ) S ) 15 S 
“25 
TS 


The following transcendental results: 


2 
S ln (28 +458 4450) = 30 
( 28<.15S 





S = 11 satiefies this expression. 
Substituting S = 11 into the expression for(pS) 
ps = 4.18 (p = 0.38) 

As a result of the above —- based on the very stringent assumption 
that search widths will vary uniformly from 0 to 30 miles (implying pools 
from 0 to 30.7 miles in diameter) - it is concluied that a search spacing 
of 11 miles maximizes the ratio of probability of detection to time of 
search effort. This corresponds to anticipating an average search width 
of 8.44 miles and a pool diameter of 9.14 miles of uniform minimm 
intensity. 

If the assumption had been that pool widths varied uniformly from 0 
to 50 miles in diameter the conclusion would have been that a spacing of 
25 miles is ideal. It is apparent that the conclusion is very sensitive 
to this assumption. Statistical data are required to provide a firmer 
base for this conclusion, but until data becomes available the assumption 
as stated is considered reasonable. 

b. Small location Error. 

When the standard deviation of the location error is only 5.6 
miles the effect of the current is relatively significant. After a delay 
of 24 hours the peak of the distribution would have moved out to a circle 
of 12 miles radius (assuming a drift of 0.5 knots). The search path is 
along the peak at the time of arrival (modifying spiral to an expanding 
square). After the first path the,"a posteriori" probability becomes 
bimodal with a circular peak inside and one outside. The pattern con- 
tinues to search these maxima. 

It is beyond the scope of this thesis to develop a probability of 
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detection statement for this situation, but it is worth considering the 
tendency of the "a posteriori" distribution. As each peak is searched 
the surface of probability becomes flatter until after one pass of the 
area the distribution approaches a uniform one. As the first search is 
being made, however, the situation is not much different from that de- 
veloped in the case of the larger deviation. In the latter case the 
probability peak was continuously being pushed outward. In the situation 
of the moving target the circular peak is being split into two concentric 
rings enclosing the original peak. Each resulting ring is searched, and 
the peak is pushed both inward and outward as the search progresses. 

In order to get some feel for the problem (which is more difficult 
mathematically than physically) the same expression for the probability 
of detection is used as a foundation for investigating the search width. 

From a geometry standpoint, after an area of a given peripheral 
size is searched, 0.75 would receive the same intensity of search in the 


current problem as was received in the stationary, large deviation prob- 





lem. The remaining quarter would be only half covered but with twice 
(i a ee - 





the effort. 

: 
STATIONARY TARGET, MOVING TARGES ————---— 
LARGE DEVIATION SMALL DEVIATION 
P= ~(1+nW exp (pW ) F20.75) 1- “Cag exp td 40.125) 1 1-(4ze) 

me) $s") S) 
exp(-2W ) 

Be 1 (Ss zm) 


Fig. 5 Comparison of search patterns. 
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Under these very rough approximations and the assimptions of the 


previous discussion a table similar to Table 2 can be computed. 


Associated Spacing Average Prob. of Measure of 
Pool Diameter Cy. Detection (p) Effectiveness (BS) 
2 2202 -767 1.55 
4 Rees °70 2025 
6 4.08 0654 2267 
§ 4-19 soto 2095 
10 5 041 582 315 
20 Vel? e472 3.68 
30 9.61 e409 aes 


TABIE 3. Comparison of spacings using the MEASURE OF EFFECTIVENESS (SMALL 
LOCATION ERROR). 

This indicates that the largest spacing would maximize the appraxi- 
mate probability of detection per time of search effort provided the 
actual probability of detection is at least proportional to the above 
estimation and that pool sizes are equally likely in size from 0 to 30 
miles. 

After one pass of the area using the doubly expanding square and a 
spacing of 9.61 miles the target location distribution approaches a UNI- 
FORM one. It is obviously not UNIFORM, but the peaks have been flattened 
and multiplied and the continued movement of the current throughout the 
search builds up the probability of the target being in a path already 
swept. The end result is that the unit mass is best considered to be 
evenly spread over the area. 

As a result of this reasoning a parallel sweep search will maximize 
the probability of detection per time of search effort after the first 
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pass. The expression for probability of detection then becomes: 


P=W for SW; P=1 for SW, 
S 


Following the same approach as in subparagraph 11.a. it is desired 
to determine the spacing that meets the established measure of effective— 
ness - again under the assumption of equally likely search widths (pool 
diameters ) over the range of 0 to 30 miles. 


30 
Dp = st 2 dw for SW 


slit 


We are interested in pS. Since p cannot exceed 1 there is no sense 
in making S smaller than 15 because p would remain constant as S was made 
smaller and pS would therefore decrease. Obviously the maximum is reached 
when S = W. Therefore this search should be conducted at a spacing of 15 
miles. 

For the second pass the search area should be a square whose dimen- 
sions are dictated by: 

a. Three times the location standard deviation to provide an almost 
certainty that the pool was included in this circle at the time of the 
fix, plus, 

‘be. the product of the drift and the time that would elapse between 
the time of the fix and the time of completion of sweep of the area on 
this second pass. 

The sweeps should start along the side perpendicular to the most 


frequent direction of the current at 15 mile spacing. 
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12. Aircraft capabilities. 

It is assumed that an aircraft will be available for the search 
equipped with the proper survey meter and water sample floats capable of 
the following: 

a. Maximum range to scene 1000 miles with minimum of 3 hours on 
station available. 

be Search speed 180 knots. 

13. Search procedure. 

This problem deals with an underwater burst which vents to the ate 
mosphere. All of the effort of the study has been devoted to an analysis 
of activity at the scene, but the study has presupposed that an air 
sampling barrier would be established downwind from the blast. Any 
samples obtained from this endeavor would not be subject to masking am 
quite informative about the circumstances of the origin of the blast. 

As previously mentioned, correlation of hydroacoustic information 
may take several hours. It is reasonable to assume that as much as 24 
hours may pass before search craft could arrive at the scene. The appro- 
priate search plan should be followed depending on the physical situation. 
If a pool or suspected pool is located an attempt should be made to locate 
the area of highest intensity (with due regard for safety of the crew). 
Intensity readings should be taken and correlated with position and time 
as accurate as possible. Relative positions within the pool are more 
important than exact geographical locations. At the point of highest 
intensity water sampling floats should be dropped for later recovery by 
surface craft or helicopters. 

The search aircraft should remain in the area as long as safety and 
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operational requirements permit so it can direct other searching activity 
for recovering water sample floats. 

Frequent intensity readings may provide sufficient information to 
establish decay curves to extrapolate back to the time of detonation. 
Great accuracy would be necessary in correlating the movement of the pool, 
height of aircraft, location of aircraft over some relative point of the 
pool, etc. to permit much of an estimate of the decay rate. The best 
procedure would be to keep recording the intensity of the area of highest 
intensity to provide some semblance of a reference. This would corrobo- 
rate hydroacoustic information. 

Normal search procedures should be followed in marking datum and 
search area, and making reports. 

14 Conclusions. 

The problem was to determine the optimal searching procedure for 
locating the site of an underwater explosion which had vented, in order 
to classify it as nuclear or non-nuclear. 

A hydroacoustic network is anticipated to be the detecting agent. 
Such a neiweck has no ability to distinguish an underwater nuclear explo- 
sion from a non-nuclear one. Classification must be accomplished by a 
thorough investigation of the scene for evidence of nuclear activity. If 
no such evidence is found a probability statement can be made of the 
nature of the explosion, and upon this information a decision to terminate 
or continue the search can be'reached. 

A MEASURE OF EFFECTIVENESS was established to be: 

Maximize the probability of detection per time of search effort. 

Several parameters were investigated to appreciate their effect on 
the problem. 
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a. Height. 
It was determined that the maximum intensity reading at the 


meter occurs on the surface. As altitude is increased the intensity 
decreases up to an altitude above 20,000 feet. After that the increase 
is gradual and almost negligible. Of course at such a height the in- 
terest is only academic since very few waves have sufficient energy to 
survive that thickness of air. 

b. Search Spacing. 

The development for search spacing as presented by Koopman /2/ 
was expanded and adopted for this study. 

ce. Search Width. 

The effective search width was determined to depend on the radio- 
active pool diameter which was an unknown quantity. Pool diameters depend 
on yield, depth of detonation, and physical characteristics of the water 
and the exploded device. In order to discuss the problem quantitatively 
an averaging procedure was used, assuming any pool diameter to be equally 
likely between O and 30 miles. This is a very brash assumption, and the 
results of the study depend heavily on this weighting function. Unfortu~- 
nately, little statistical data are available which might indicate the 
distribution of this function. This assumed UNIFORM distribution of 
weights providesgan unexpected conclusion when the criterion is the MEASURE 
OF EFFECTIVENESS mentioned above. If the weight function favored some 
particular pool diameter it would be reasonable to adopt the most effect- 
ive spacing for that diameter; however if there is no bias toward a par- 
ticular width the conclusion is not so evident. Various spacings were 
concluded to be optimum depending on the standard deviation of the dis- 
tributions considered. 
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d. location Error. 

There are two categories of location error considered. large 
errors with standard deviations above 25 miles dictated a search plan 
that ignored current drift. For smaller deviations the drift had to be 
reckoned with in designing the search plan. The location error consid- 
ered the effect of both the initial hydroacoustic fix error and naviga~ 
tional errors of aircraft arriving at a given spot in the open ocean. 

When current movement was considered (with small standard deviations 
of location) the direction was assumed to be random, but the drift can be 
predicted within an acceptable margin. 

e. Search Plans. 

When the standard deviation of the location distribution is 
large the squares of uniform coverage maximize the probability of detec- 
tion per time of search effort. On the other hand if the deviation is 
small,a retiring search maximizes this criterion for the first pass. For 
subsequent passes of the small deviation problem the "a posteriori" 
rippled distribution approaches a UNIFORM, and a parallel sweep search is 
considered best. 

Finally, although this comment is beyond the scope of the stated 
problem, it seems worth mentioning in passing that very little special 
equipment is required to perform this function of classifying detected 
disturbances. Water sample floats which can be dropped from an aircraft, 
sample the water and seal themselves — to be recovered by surface craft 
or helicopters, and the survey meters described by NRDL /1/ are the only 
non-standard items needed. These items can be based ashore and on 
carriers and need not be on aircraft until a plane is assigned a particu- 
lar task of classifying. 
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APPENDIX I 
HEIGHT OF SURVEYING AIRCRAFT 

In order to determine an appropriate search plan the search height 
must be determined. According to reference /1/ the meter will be con- 
structed with a shield forming a cone about the axis to prevent excessive 
interference from cosmic radiation. This shielding forms a central angle 
of 150°. 

As the search craft increases its height the area seen by the meter 
will inerease and one would expect this to increase the intensity at the 
instrument. However, the attenuation due to the atmosphere will decrease 
the intensity at the meter. 

Model assumptions: 

1. The pool is UNIFORM and INFINITE in diameter. 

2. Gamma énergy at the surface is 1 MEV and the "build-up" 
factor is not considered. 

3. Absorption is linear varying at a constant rate with dis- 
tance from the source. 

4. Height of meter is much greater than the diameter of the 
meter opening. 

I = Intensity of radiation at the meter. 

A = Area of meter opening. 

I = Intensity of point sources at the surface. 

h = Height of instrument. 

u = linear coefficient of absorption in air (2.265 x 10>/FT), 

x = Distance of meter from radiating source. 

r = Distance of source along water from point on surface directly 

below aircraft. 
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Angle of rotation of "r" on surface. 


g 


w= Solid angle subtended by meter opening at radiating source. 


Angle that point radiating source makes with meter. 


The portion of the radiation that arrives at the meter from each 
point source is equal to the ratio of the solid angle (w) subtended 


divided by 4. 


W A 
pom 


x 
This assumes that h is much greater than the diameter of A; there- 
fore the projection of the area of the meter opening on a plane tangent 


to a sphere expanded about the point source is approximately equal to 


Peajectean 
© ee ov Qe o 
\Low sphere 


the projection onto the sphere. 





Fig. 7 Parameters of integration. 


htan @ ;2// 
i= I A cos g exp[ —ux | r dr dO 
Jf ane! 
0 O 
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(x2-n2 )e 


cos J =h r= 
x 
dr = x 
x<<=h 
h sec 9 Ls 
a) 61 -f Te A b spp{ uy dx 
h 


It would be of interest to find a height which would maximize the 


intensity at the meter if such exists. 
The integrand in (1) is continuous except at o and the derivative 


can be investigated by applying Leibnitz's Rule: 


(2) aL = ath a B) sec P - b exp{mu h) 
f sec ff es -} 
+ sap \—wx dx | 
~~ 


Consider the last integral in (2). 


h sec # “Th sec & 
a spor ky sanp( na sal | 
s x x x p 


h sec @ 
= | -exp(-ux eu (In xauxtu223x7 tute, 
x 2°21 3°31 AeA! 
This term by term integration is justified by virtue of uniform 
convergence of the infinite sertes expansion. 
Substituting this in (2): 
(3) aI = -In sec # + u h(sec $-1) - uth? (sec* f-1) 
dh 4 
+ wh? (sec? 1 ) = yAn4 (sec4—1 ) ee 
18 96 


Sec @ = 3.8637 


oo 


dI was programmed on the 1604 computer for a large number of values 
dh 


of h. The results showed that the derivative was always negative, be- 
coming smaller in absolute value as height increased. This indicates 
that a plot of the intensity would be an exponential decay with its maxi- 
mum value at h = 0, 

Assumption 3 in the previous model is not realistic because the 
attenuation varies with density. Since air density varies with height 
the coefficient of absorption must also vary with height. The mass co- 
efficient of absorption is equal to the linear coefficient divided by 
the density and it does not vary with the chemical or physical condition 
of the medium. For air this quantity is 0.0274 om*/gm. 

According to the Barometer equation derived in STATISTICAL 


MECHANICS 


P = Py exp(-mgh) 
KT 


oS 
ll 


Density at height of instrument. 


m = Mass of a molecule of the medium. 


g = Acceleration due to gravity. 

h = Height of the instrument. 

K = Boltzman's constant. 

T = Temperature (degrees KELVIN). 

P,= Atmospheric density at sea level. 


u'= Mass coefficient of absorption. 


mo 
i 


Slant range from radiating source to meter. 
§ = Limit of shielding (75°). 


Returning to equation (1) and substituting: 
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exp (-ux] = exp|—u" Px| = exp|-u'xP, exp(-aah ) | 
h sec 9 
a le 5 h exp|-u'x?, exp(-mgh) [ax 


pleg'P, exp(-gh) h 'P exp(-mgh)h 
= sxpl- © exp(-meh) sec 6 | sxo[-u' Py exp(-aeh)h | 
a 


(4) 


wi- 
ae 


sec 9 
h mg exp(~mgh )u' Pa exp| -u'P, exp (-mgh) x|ax 
KT KT x KT 


h sec @ . 
+ expj-u' Py exp ~ngh) x fax B = IgA 
f 2. KT 
h 
The last integral of equation (4) becomes: 


h sec 9 h sec 
(5) |= exp(~ax ) € exp(-ex) dx 
x x Where a=u'P, exp (—mgh ) 
h h KT 


The first part of (5) balances the first two terms of (4). 


h sec @ ° 
a(h mg -1) f exp (-ax )dx 
KT 
h 


1 ai 


B dh x 


(6) 


a(h mg <1 )| a sec § - ah(sec f-1) + ach? (sec~ f-1 ee 
KT 2°2) 


Consider the following: 
"a" 4s positive for all values of h; 
The infinite series is always positive; 


hmg=-1 is O when h = KT? 
KT mg 


Therefore there is a maximum or minimum or flex at h = KT 
mg 


(= 3d sec Pea h(sec Z-1) + ah? (seo J~1).. i 
202! 


a°I 


adh2 th = 





a 
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Again "a*, B and the series are positive. The conclusion is that 


h = KT is a minimm point. 
mg 


Equation (6) was programmed out to 30 terms for values of height 
from -100 to 40,000 feet and several sensitivity checks were made. 

The temperature is recognized as a variable over the range investi- 
gated and it was varied from 300° (81°F) to 273°K (32°F). The values 
varied slightly but the structure of the curve remained the same. A few 
points were checked with the temperature at 243°K (-40°F) and again the 
height of the minimum decreased slightly but the structure of the curve 
remained the same. This covered the range that average temperatures 
could reach and the minimm varied from 21100 feet to 26070 feet. 

Variations in surface air density (Ee) were looked at but their 
effect were even less significant. 

Variations in "g" were not looked at, but they are obviously of no 
concern. 

Lastly, the evaluation was made for 10 terms of the series and then 
for 30 terms of the series. In both cases the values were the same out 
to eight significant figures. 

As a final investigation the expression for the intensity was pro- 
grammed and curves obtained are presented in figure 8. It is satisfying 
that the same MINIMUM points were reached by two independent approaches 
and two different programs. 

The development was as follows: 
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h sec @ h sec @ 


= Igh| -exp(-ax) “4 exp(-ax) dx 
2 X x 
h h 
= exp(-h sec 9 = gga) a sec fa h(sec P-1)... 
2 jh sec 


LoA (1 Lira fan sec $+ (1-1 1_)ha (sec J~1 ) 


2 ey 


{-hee 4 _4th2a-(sec® $-1)...4/+.1_- ” a (sec f-1)... 
OT eer BHT 


= Lot |.741+na(~20 sec Pr (1-4) ha(sec @-1)......) 
2 2 


In fig. 8 the bracketed quantity is plotted as f(h) - the "Intensity 
factor", 
CONCLUSIONS 

1. The greatest intensity at the meter will be experienced when the 
height is zero. 

2. The minimum point of intensity is little more than academic 
interest since the assumption of a UNIFORM source of diameter at least 
equal to the field of the meter becomes a stringent condition at 25,000 
feet. The pool would have to be 30 miles in diameter. More directly, 
this problem is interested in maximum intensities. 

3. Operational conditions will dictate the height of search with 
the aim of flying as low as feasible on a prolonged mission. Figure 8 
demonstrates the degradation of intensity at the meter (and therefore 
the decrease in search width) for the height of search selected. It is 
not possible to predict the decrease in probability of detection which 
will be experienced as height increases because the physical conditions 
of the pool of radiation are not know. Since the meter has a threshold 


intensity which must be reached before a signal can be recognized, flying 


Be) 





ata greater altitude than operationally necessary may Aégotarcta ti 
intensity at the meter so that the threshold is not reached - whereas 
it may have been reached at the lower level. 

4. At the altitudes that are operationally feasible the curve will 
give a reliable indication of Intensity factor (f(h)) for any atmospherié« 
conditions (neglecting the build-up factor effects) since it is insensi- 


tive to temperature changes below 3000 feet. 
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APPENDIX IT] 
SEARCH THEORY 
This appendix is based on Koopman's Search and Screening (OEG Report 
No. 56) /2/. An attempt is made to present his development in more de- 
tail as it applies to this thesis’,. 
A. Development of probabi Wie of detection of a target of lateral 


range X from a straight aircraft track: 
(1) p(X) =1- exp[-0.092(z)? | 
x 


1. Consider the case of "continuous looking" where "s" is the 
instantaneous probability density of detecting a target. If s remains 
constant throughout a search p(t) = 1-exp(-st) where p(t) is the prob- 
ability of detection during time t. 

If s is a changing density because of relative movement with respect 
to the observer, and is time - dependent, the expression is more compli- 
cated. The instantaneous density s is inherently dependent on range and 


the fact that it varies with time is conveyed by the symbol s,. 


p(t) =1 - oxp|- [x at | 
O 





Cbser nec 


Fig. 9 Target's relative track. 


: 
The target is at (x,y) at time t and x = x(t), y = y(t). Let the 
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initial position at t = t' be X, = X(t"), y, = y(t"), and the final posi- 


tion at t = t" is X = X(t"), y, = y(t") with the target tracing a track 


Wl 


along C in figure 9 above. 


# = o| (x2(t )+y2(t) 8 | = St 


The probability of detection along path © is: 


+" 
Pe =i- ex I a ((e2teey2(e Pat | 


The integral is a line integral along ©; if "w" is the relative 


speed and "1" is the arc length along © from X,, Y,, then 
Eo | exp|- fs ({x2(t)+y2(t J] 4) dl aige=ada: 
c WwW WwW 


Following Koopman's notation the exponent is called the sighting 
potential F{C]. This quantity possesses the property of additivity so 
that where C = C, + Co, then F{C] = F{C,} + Fic] 

In the case where X is the lateral range, X is constant and y = wt 


(w is a constant), 


+" y" 
F)c} = f s( (x24w2t2|F)at =1 J 8 ( fx@4y7]? day 
W 
t! y! 
According to the inverse cube law of sighting where h << 


(2) s=kh 
r3 


{ql - it 
a §? J (F497) 3/2 
y' 


‘ 
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Fig. 10 Detection at fixed speed and course. 


When the sweep is "relatively" infinite 96' = 9" = 90° and - 


3) F(CJ=FLxy= ae Gri)s enh 
—— 


When the observer and the target are on straight courses at constant 
speeds for a long time before and after their closest point of approach, 
the probability of detection is a function of the lateral range X. 

fh = [~ £¥° [- Fc x 

2. Parallel sweeps - 

A target is presumed to be at rest on the ocean in an unknown 
position, all equal areas having the same chance of containing it. A 
search is made along a large ("infinite") number of parallel lines a 

AL 





common distance apart. This distance is referred to as Sweep Spacing (s) 
and the problem is to determine the probability of detection of a target, 
P(s), by this configuration. Other questions concerning probability of 
detection, if target speed and course are known or have a distribution 


of their own, can follow in development. 


Consider the lateral ranges of the target from various sweeps 


(x, 4) 


a near wwe meses 


ao = © = Pamee 
Fig. 11 Parallel sweeps. 





For paths to the left the distances are X, X + 8s, X + 2s... and for 


those to the right s - X, 28 - X, 3s - X... Both cases can be combined 


in the absolute expression 
lateral range = |x - nf | where OC XS S§ 


n= 0, £1, £2, £3, one 
The probability of detection by the n*® lateral sweep when the tar- 
get is at (X, y) is: 


Fy (x)= j- ere L- FCI K~ ns) | 


The probability of no detection by the nth sweep is 1 - P,3 and no 
detection by any sweep is the infinite product TT (1 - p,) for all values 


of n. The probability that at least one sweep detects a target at 


ory) is 


i) ple slate Texel F(vemills - exp pis] 


— 
— 


ae 4, 
(y) where Pixs) = = FAYx-ns\) ; Os “CS 


/ haeee 
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Consider the lateral distance of the target X as having a uniform 
distribution between 0 and s. P(s) is the average of P(X, s) over all 
values of X in the interval. The expected or average value of the 
probability of detection for a given s can be expressed as a function 


of s. 


3 ; cm ov 
P? Pos) PC Ss dx = z Sa- exe [- dixs) Jelx 


This is the solution to the problem but it remains to express it 
more quantitatively. 
The effective visibility (E) is defined as half that sweep spacing 
for which the probability of detection by parallel sweeps is 4; i.e. 
P(2E) = 4. 
‘In the definite range law detection occurs if and only if, the tar- 
get happens to be within the definite range R of either of two adjacent 


sweeps. The chance for this is 2R = W when s >2R = W, and unity when 
S 9 


s<S W. From this it follows that E = W (-since W = .5; W. = .5; W=E), 
8 25 


For the inverse cube law when h << r we get from equation (3): 


F(X) = 2 hk; and from equation (5): O(X, 8) = 2hk 
xe a 


This can be restated as # (X, s) = = 7 ese“7X because of the 
g S 


following development: 
a. The infinite oe oo ae for sin 4 is 
sing= 2 T (I~ y 
on. 
b. Take In of el side - 
Insing= 4, Zt 24(\- S74) 
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c. Differentiate both ao - 


cos 2 _ 5 Z 45 
sind . : ih a we 


ad. Which may be rewritten: 
Cx? 
| = meee 


cot 4=-=—=- ~ 


= 
x fd- (14 “Tl qt tee zh) 
Note: This is the Fourier expression for cot 4 and could have been used 
as the starting point of this development. 


e. Break up d. and worre for A and B by partial fractions: 


cot 3 = a 2 fee alse ) 
.* qe\\ wire v varie / 
Bs 
From partial fractions = 
A(nll+ 2) + BlnTT- 3) - 22 
AS - B& = 26 
A=2+B 
An T+ Bnff= 0 
2+B+B=0 gore B= =1 
A = 2 +9Baa4 
a 
(7) ot 4 “me = 2 (+> h T- 42 NT + % 


te-+- on 
() ct 2@-i+ S 
tc ne-0? Beni 
2325 hx oO 
Orscole= 2 ales 
how Om ge rir 
The only way to see the above steps is to substitute values for n in 
(7) and see that you get two terms for each n which can be represented by 
the expanded limits of the expression in (8). Note also that when n = 0 


in (8) the missing term is 1 which has thoughtfully been provided. 
- 


f. Differentiate (9): 
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one 
a... 
oe Coot =. (Z-nr)* 


g Let & = Te and weir Teee (10) 


1 ee = + 
Cb Ga 
(41) Then (h(X%,S) = hk & ~ 245R 7? IT 
(p S) See a ae Lis oh eas ae 


Inserting this in (6): 
-_ {° RA Te ech EX) TS 
P(s) = sf [i-exP(~ 22% We ese* 2) Jon 


There is one more a of wizardry which Koopman gd 


Ur S ‘ 


ate) = PCs)= 2 Pada ra Css “eg / 
Differentiate with respect to B 
n “- 


O : | 
4 eye (-8) CxP(—Beot*s) cse*s do 





I} 


ese°@ dO = -d(cot 9) 


when © = 0 cot 0 


0= TT cot @ 
2 


I 
© 


Let eot~o = x? 


d cot 9 = dx 
re 


g'(B) = — + @ xP (-B) [exeCe mad 
& exp(-8) + L (Tee exe(-B) 
° (er) 


Integrate #'(B) with respect to (B) from 0 to B! (a particular value 


of B) 


48 





Let X* 


= B 2 X dx = dB 
xX =Be when B= 0 X = 0; 
2 when B = B!, XxX = (P yz 
Q(B) = exe (-8) 2 (er (~ x3) 2yad 
J an 78 - Ga “SP : 


x | 
(12) O(B) = An f ey 2l~ ie) dy 
(pl 2 = x! = a particular value of X. 
By definition (12) is the error function (erf) based on the normal 


distribution. 
Hence P(S) = O(B) = erf (Bye = erf T fan J? 


8\wW 
For the inverse cube law the effective sweep width W = 22 ka)? 
Ww 
(13) P(s) = erf W(T * 
8 2 
P(2E) = 4 = erf (Tye w 
LE 
Since erf 0.477 = 0.5, then (JEW = .477 
LE 
and W = 1.076E, hk = .O46E* 
Ww 


Substituting this in (13) 


P(s) = erf (0.954E ) 
8 


It is of particular interest at this time to note that the proba- 
bility of detection is not particularly sensitive to the law by which 


the spacing (s) is determined ~ particularly at lower probabilities. 
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Fig. 12 Probabilities with parallel sweeps. 
Recall the following equations: 
p(x) e-1< exp[-F (x ) | for lateral range distribution ; 


fad, ) p(X) =1- exp “2s, for inverse cube law ; 


W = 1.0765 from the error function development - 


(15) hk = 0.046E* from the error fumetion development 
ro 


Substituting (15) into (14): 
P(X) =1- -0.092/E\" 
(X) exp] - 9 (BF 


Which is equation (1), the immediate object of the above expansions. 

Before proceeding, let us review the principal assumptions upon 
which the "error function" development is based. 

a. Target at rest in an unknown position with lateral distribution 
UNIFORM between sweeps. 

b. Observer is on a straight course at a constant speed for a "long" 
time before and after the closest point of approach to the target. (Sweep 
is "relatively infinite"). 

c. Inverse cube law of sighting applies and hvcr. (h" is the 
vertical height of the observer and "r" is the horizontal distance from 
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observer to center of the target. ) 
d. Sweeps are independent. 


B. SPACING (S), 


me Consider the SQUARE OF UNIFORM COVERAGE search pattern (section 


7.3.2 of reference 2). 





OWT Re AL 


Fig. 13 Square of Uniform Coverage. 
let the normal density function o 
| oe 
(16) COs a)=F Cr) = spee OW OTe 


represent the probability that the target at the time of the fix be in 

the small region dx dy at the point (x,y), a distance of r from the origin 
and where 6 is the standard deviation of fix accuracy. The first problem 
is to find a spacing so that on the initial square the probability of 
detection per unit time shall be a maximm. 

As was indicated before, the spacing is not sensitive to the law of 
detection, and the inverse cube law development will be followed so that 
the probability of detection of a target of lateral range x from a 


straight aircraft track is 
(17) p(x) =1 - exp|-0.092(8)*/ 


x 


om 





Reflect for a moment on our problem and the assumptions of the pre- 
vious section. The same error function result is being used. 

1. The lateral distribution is certainly no longer uniform between 
sweeps; however when the standard deviation is large compared to the 
spaoing the approximation of uniformity between sweeps is not very in- 
accurate. 

2 The sweeps are not "infinite" in any one given direction, but 
except in the center the length of sweep is relatively large, and again 
the practical approximation is acceptable. 

3. The inverse cube law is assumed in the development with hc. 
First of all, as already stated, the spacing is not sensitive to the law. 
Caution must be used in proceeding, however, because although the distance 
to the center of the detected configuration can normally be expected to 
be greater than the height of the aircraft, in the case of an extremely 
weak source approaching the minimum detection capability of the meter = 
when directly over the pool - the range to the center would be less than 
the height and the probability of detection would be overestimated. The 


instantaneovs probability (s) would equal kh vice and the error 
Gzng3/2" °° 3 


would increase as the ratio of height to horizontal distance increased. 
4. Sweeps are still considered independent. 
The probability of the target lying on the strip parallel to the y 
axis between x and x + dx is found by integrating equation (16) over all 


values of y and obtaining the marginal: 
J (Fe yay s eon exe -S,) Ay 
qn K din = je So (~ 4° ) 
(18) AK) Teng o( 211) a¢* 
The peak of the distribution is located at the center. If an 
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indefinite straight flight is made along the y axis and has failed to 
detect the target the distribution is altered and the differential co- 
efficient of (18) no longer represents the lateral density of the tar- 
gets. By Bayes theorem the "a posteriori" probability considers the "a 
priori" probability of (18) above and the "productive" probability of 


not detecting the target on a sweep through the center: 
2, 
l— PCY) x OxXe {- 0.052 es. I (See equation (17)) 


Bayes Theorem states: 


Let the union of events B, be the sample space where Bs 


a 
(4 = 1, 2, ..., k) are mutually exclusive and exhaustive. 


Let event C also be defined on the same sample space 
P(Ba\C)=  PCCIBa) PC Ba) 
& P(clBr) PCB.) 
Ca 
In our problem let: 
GC = not detecting target 
By= target located in strip (x, x+dx) 
The denominator of Bayes theorem can be expressed in the continuous 
case as Oo 
f Cow PCC/4) C4 ) A rai where y represents any strip, 
This integral is just a mumber, say K, which actually normalizes the 
numerator so that instead of speaking of the probability of the event Bs 


given event C we speak of the density function 


fin): 3 (8 41¢) e = ribs) 


- so that the new density function is proportional to P(C |B, )e(B,) and 


this translates to 
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(19) Liq Peas £09) = ki ex’ /-, - p78 (5) 
a ‘WS 


in words - the probability of not detecting the target given that it 
is in a strip (x, xtdx) times the probability that it is in a strip 
(x, x+dx ) is proportional to the new density of the probability that the 
target is in a strip (x, x+tdx) given that it was not located on the first 
sweep. 

To find the maximum of the new distribution differentiate (19) and 


solve for the x which makes the expression zero. 


? x. oe 
—. -~ 260.094) EF 
ae ee 


x? 
w= | Oo. oeeE CE s)'* 

So the maximum probability target position is no longer humped at 
zero. Double peaks now occur with a depression at the origin. However 
the distribution is skewed, and the distance between the first and second 
tracks, D, should be slightly greater than 65(EN2 to cover a larger 
volume of probability on the second sweep, 

The probability of not detecting on the second sweep given that the 
target is ina strip(x, xtdx) times the probability that it is in a strip 
(x, xtdx) is proportional to the new density of the probability that the 
target is in a strip (x, x+dx) given that it was not located on the second 


sweep. This follows from the same reasoning used in equation (19) 


te 
ECan) Ax = exp l- ony ~ 0.042 LE am a2 ie E* [ay 
| / ‘alin (X- O) = 
Integrating this function for all x and ‘solving for a D which would 
minimize the remaining function we would be maximizing the probability 
of detection on the second pass. Koopman presents the solution, arrived 


at by numerical integration, in section 7.3.2 of SEARCH AND SCREENING /2/ 
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as 0.75 (ES)2. This spacing is then used for searching about a point of 
Pix. 
C. Search Patterns 

1. Square search for a stationary target (Square of Uniform Cover~ 
age"), 

The details of the theory are described very explicitly in 
SEARCH AND SCREENING /2/. This search assumes a circular distribution 
of target position as a result of a probablistic fix. This pattern 
attempts to give a high probability of detection per unit time at the 
beginning of the search when covering the area where the target is most 
likely to be. 

2. Retiring square search for a moving target (Section 7.3.3 of 
SEARCH AND SCREENING /2/). 

In this development it is assumed that the direction of movement 
4s not known but uniformly distributed (1.e., equally likely in any direc- 
tion), while the speed of movement is reasonably well established. 

The theory is to search the peak of the distribution at the time 
of arrival. It has moved outward from the fix a radius of Ut where U is 
the speed of target movement and t is time late. During the search the 
target continues to move and the pattern is to search the peak, then out- 
side, then inside, then outside, ... etc. This will give the maximm 
probability of detection per unit time of search effort. 

It is recognized that both of these patterns present operational 
problems in implementing. Expanding searches are always frowned upon 
because of the navigational difficulty in making legs precise. This is 


particularly true when spacing is small - less than 5 miles. 
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APPENDIX III 
SEARCH WIDTH 
USNRDL Technical Report 525 /1/ envisionsan airborne detection in- 
strument capable of sensing radioactive surface activity on the sea when 
the intensity at the meter is 0.2 microroentgens/hr. At a search height 
of 500 feet the intensity of a broad source at the surface would have to 
be 0.74 microroentgens per hour. The meter envisioned would be colli- 
mated about a central angle of 150° to reduce cosmic interference. 
A. For determination of the effective search width consider the fol- 
lowing model: 
1. Physical parameters: 

a. Search height — 500 feet. 

b. Meter described above (1.¢e., shielded around are of 150°, 
and capable of detecting 0.2 micror/hr). 

c. Minimum surface intensity present - 0.74 micror/hr. 

d. Meter system time constant about 2 seconds. (The technical 
report anticipates 0.2 sec for the meter and it is presumed that a re- 
corder, alarm and other auxiliary equipment would make system time con- 
stant an order of magnitude greater). 

2. Assumptions: 

a. Surface pool of radiation is circular and of uniform inten- 
sity. 

be. The radiation on the surface ends abruptly at a definite 
boundary. 

c. Definite Range Law of detection obtains (W = E). 

Due to the geometry the meter will see a circular area about 0.613 


miles in diameter (See Fig. 14). 
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Fig. 14 Meter geometry (side view). 

If a circular area of 0.74 micro-r/hr intensity and 0.613 miles 
diameter were present in the search area, the effective visibility (E) 
would be zero since the plane would have to fly directly over the center 
in order to detect it. 

If the pool were larger, the search width could be considered as the 


diameter of the pool minus the meter visibility (0.613 miles). 
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Fig. 15 Surface geometry. 
Another consideration in determining the effective visibility (E) is 
the meter system response time. In the USNRDL report /1/ the following 


standard formula is provided. 


Di 





1) 626, - exp (+t /T] ) 


where S = meter indication 


So= source strength 


t 


It 


time over source 
T = system time constant 


Neglecting attenuation and build-up, if the meter were over the 


source for ten seconds the meter would indicate over 99% of the source 


strength. 


For an aircraft traveling at 180 knots the meter would have to be 


over the source for 0.5 miles. This cuts down the effective diameter 


another 0.1 miles for a pool of two miles diameter and a lesser amount 


for larger pools. 
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Fig. 16 System time constant effect. 
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For this problem a correction of -0.7 miles is made to pool diameter 
to obtain effective visibility. 

Since the size of the pool is an unknown in the search several sizes 
will be looked at to determine sensitivity. 

For a search of a stationary object using a Square of Uniform Cover- 
age pattern assume a fix error (2 standard deviations) to be 50 miles and 
an error in navigation of 10 miles in locating point of fix. The total 
deviation would be 25.5 miles. For simplicity we will adopt the follow- 
ing expression for the probability of detection: 


ie) p = "= (14g) exe “at (See Koopman /2/ pg 113 and 114) 


n = # passes over area; W = search width = E; S = spacing. 


ob. of detection. 
















W 
S 
2 1.3 | 4.31 ; «302 
4, 3.3 6.86 | 481 
é 5.3 8.7 | .609 
g | Wed 1Oee ) 2715 
10 | 9.3 T1135 | 809 193 
20 | 19.3 16.6 1.16 i 6321 | .673 | 857 
30} 29.3 20.4 1.435)| 42 | .752 | .928 
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TABLE 4. Uniform square coverage probabilities. 

Since it is unlikely that data are available for such small intensities 
(0.74 micro rfnr), it is desirable to consider intensities that are more 
likely to be recorded and scale effects from available search data. 

Consider a circular pool of higher intensity, say 0.05 r/hr. If the 


attenuating distance of two configurations were equal, then for each to 


59 





induce the same intensity at the meter, the product of each area and sur- 


face intensity must equal the same product of the others. 


6) 1,4,-1,A, I, = Minimm detectable surface 
6, Te x10" rT (0.3065 )= oS A, intensity when spread 
f= 4-87T WI” sp mc uniformly over A, 
I, = Uniform intensity of new 
configuration 
A, = Area of new configuration 
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Fig. 17 Geometry of detector view of high intensities. 


Area of radiation in view of detector is the sum of the segments: 


(4) as 0° (@,-sim Gott 2 (@,- 51 &) 
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(5) 


(6) 


(7) 


(8) 


Ox 
Common chord yieldss =f — SI => 
fl 2th) a 
The effective search width = 2( facia a — ( Y; 45) 
= 2 Ch, +h) 
ee cr g ° “ek ©, 
h, | Cos 5 ) h, m is COS = 


| W=E= ZK cos ht yy Cos Se) 


From (4): 
a _ iio CO ~-~S &, = ->S 
aA at ( i Mh ) ( Oz { + 6._ ) 
From (5): A sim z =, sin es 
> 


7 
ic (\-<c0% »o )= C\— Gee @, ) 
‘he 
y= 


Equations (7) and (8) are transcendental and can be solved graph-~- 











ically. 
| Soe equation (7) From equation (8) 
miles AC 8 ACI ANS 
Pool radius || % | % |i ¢ i 
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| | om (Cont'd next page) 
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_ equation 2 (8) 
Pool radius 








TABLE 5. Data for pegnieal solution to eutivendeal Sandip 
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Fig. 18 is a plot of these curves. The points of intersection yield 
the solution for 6, 9 for various pool sizes. 

Note the limiting gase in which the pool is so large that the arc 
approaches the chord of the detector's circle. In this situation all of 
the area is described by the segment of the detector's circle. 

A = + ro* (Q - sin Qo) 

@5= 0.082 radians. 

The solutions for 92 range from .0675 radians (3°51') for + mile 
radius pool to 0.082 radians (4°42') for an infinite pool. With these 
figures and for the pools under consideration the effective search width 
varies from 

2(pool radius + .306 mi) for large pools to 

2(0.5 + .306) for the half mile radius pool. 

Therefore the search width is increased by the diameter of the de- 
tector "look" at the surface (approximately) regardless of pool size. 

Remember that a specific example of a pool of 0.05 r/hr with the 
detector at 500 feet altitude is the case under consideration. 

There is one very unrealistic assumption in the above discussion. 

The intensity at the meter is not just proportional to the area. 

The distance of the meter from the source is of principal concern ue. 
considering attenuation. Ignoring the build-up factor the interrelation-— 


ship is as follows : 
f T, Anh expl-uxJ adr de =| Le Amh eve aiadrde 
‘A 4m X? | P 41 x? 
Aun OO a of wm eye é 
(9) -2 micro r/hr = ost, Oyp (- 2¥) redrde 
> 
A, qT X 
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(The effect of the Barometer equation is considered negligible - see 
Appendix I) 

Equation (9) would be solved numerically for Ay or more meaningfully 
for the I2 required in the pool to register at the meter. 

A rough approximation for the effect of the different distances of 
the source from the meter can be obtained by comparing the intensity of a 
particle at the average distance of each configuration. 


The average distance to the particle in the case of the circle of 
minimm a is 


es : halts 0,3065 Or at mn ali he 


J CAT AS Trt 
(For circle) X = 0,223 


i 


In the case of the segment the average distance is: 

X, _ 031 78 

let I, = meter intensity for circle of minimm intensity 
Iq = meter intensity for segment 
I49 = surface intensity for circle 


. = surface intensity for segment 


Al 2X (- 444 \ vor point source. 


41 ¥° _ a 
a” f Ajo A nw by e%PC— 44 x, ) gt ®, be Tyo 0,032 Ox Cie 8¢-, o1u9) 
rope Lae Amh OP (rae ¥) TY, T,° o.ont 
A,2l2 
0.347 (0.27) <de J one [0.65 Lie 


t2¢ 
In order for the meter intensity to be the threshold intensity in 
the case of both patterns, the surface intensity of the segment would 
have to be 10.65 the intensity of the cireular configuration to compen- 
sate for the different mean distances. 
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In the previous scaling discussion the characteristics deduced would 
pertain to a surface intensity of 0.5 r/hr vice 0.05 r/hr. 

The following table will assist in scaling: 

Altitude = 500 ft; X, = .2233 u = 13.788/mi. 


ee ee 





Limiting 


Iyo(r/nr)| Ag(mi2) =| 8 ho(mi) /Xo(mi) [I29/I40 
| 40375KX107! 12.18 03064 }.3175 [10.6 & 
4.375X107° | 4.7 63063. 1.3175 {10.6 


4.375K10-5 [10.14  |.3057 [23167 [10.4 
4 o375X1074 | 21.94 305 [4314 | 9.8 
4,.375X1073 |47.73 |.2805 }.311 9.1 
4.375X1072 | 107.72 | 1809 10235} 164 








If search altitude is 1000 feet the following table can be used: 





X, = .446 
jlamiting | =|. | | | wate 
Tjo(r/nr)| Ap(mi*) =| 08 ho(mi) |Xo(mi) | 156/Tyq|Ioq «| Pood + 
1.003K1073 [10.5 |.6101 |.635 [39 95 | 1622 
4 003X107" | 22.5 6001 |.628 [34.6 173 fe 1021 
1.003K107! | 48.9 559 «1.599 419.9 00995 | 1.17 
1.003 111.2 0346 4.565 $10.45 | .0005225 0.% 





TABLE 6. Search widths for various intensities. 


66 





meter 7 


Kadvoact ive 
Pool 





MWertey Viele 


Fig. 19 Identification of symbols for previous tables. 

Interpretation of above tables: 

For scaling purposes use colum I.,, as the surface intensity which 
is known and from which it is desired to scale. 

For example, if the contour of a pool of intensity .05 r/hr is known 
and it is desired to determine the search width, enter column I59 and 
find .05 r/hr. The search width would be the pool size plus .6122 miles 
when searching at 500 feet altitude or plus 1.18 miles (interpolating ) 
for a search height of 1000 feet. 

Conclusion 

The diameter of the detector view can be added to the pool diameter 

to obtain search width when scaling from surface radiation more than one 


order of magnitude larger than the minimm circular surface intensity. 
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APPENDIX IV 
EFFECT OF OCEAN CURRENTS 
A. General. 

There are many contributing causes of currents, both external and 
internal to the water. In the deep open ocean wind is a prominent in- 
fluence. Persistent winds interact with surface water, imparting up to 
two percent of the wind velocity to water movement. The direction of 
the movement is complicated, with the resultant of the "coriolis force" 
and layer friction balancing the actual wind imparted force. 

P. H. Kuener in MARINE GEOLOGY /5/, page 29, 18 quoted in part: 


Flow in deep water far from the coast driven by a constant wind 
(1.e., drift currents) at a constand speed, the motion at the sur- 
face is directed 45° to the right (in the northern hemisphere) of 
the wind. Below the surface the direction gradually diverges 
farther and farther to the right of the wind and concomitantly the 
velocity decreases. At a depth called the friction depth, which 
varies with the latitude and force of the wind and does not exceed 
200 m, the current is directly opposite to the surface current and 
has fallen to 4% of the surface velocity. If bottom friction is 
left out of account the average movement is at right angles to the 
direction of the wind. The frictional force exerted by the wind 
and the coriolis force at right angles to the movement are then 
equal and exactly opposed in direction. This state is closely 
attained in nature. Yet there is some slight friction, and the 
resulting force plus the "coriolis force" is equal and opposite to 
the force exerted by the wind...... 


When the wind just begins to blow the first movement of the water 
will be in the same direction as the wind, but as soon as motion 

of the water sets in the rotation of the earth makes itself felt. 
In shallow water, where the friction with the bottom is consider- 
able, the current will be less deflected by the coriolis force than 
in the ideal ¢gase. In regions like the North Sea and the Baltic, 
where the winds are variable, where the water is shallow, and where 
the surrounding coasts are at no great distance wind driven currents 
will generally follow the direction of the motive force fairly 
closely. 


Internal characteristics and variants create an isobaric gradieht 
between waters. Salinity, temperature, pressure, level, etc. all have a 


bearing on currents; however, considering all open ocean currents, none 
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are stronger than the GULF STREAM which has a velocity of up to 3 knots. 
Most open water movements run between 0.3 and 0.7 knots, but it is very 
difficult to predict the direction of flow at any given location at any 
given time. 
Bowditch /4/ states: 
Current arrows on nautical charts represent average conditions and 
should not be considered reliable predictions of the conditions to 
be encountered at any given time. 
B. Analysis of Currents in Sample Ocean Areas. 
1. Set. 

Enclosed are data on current direction and speed. The first 
eight (8) pages (figures 20 - 27) are plots of the frequency of current 
directions (set) for the specified areas for the months indicated. Data 
were collected by the Hydrographic Office from ship reports. The number 
of observations upon which the graphs are based is indicated inside each 
curve. This figure is followed by a single number representing the per- 
cent of observations which reported no current (no set or drift). The 
area represented, broken into 5° X 5° quadrangles, is a ten degree strip 
of latitude (25°N to 35°N) across the Pacific from 150°E to 170°W longi- 
tude (about 1+ million square miles). Choice of the area was influenced 
by the large number of observations available over such a large expanse 
of open ocean. The months plotted were chosen arbitrarily to provide a 
three month period to observe gradual changes (January, February, and 
March) and other samples during the year to show overall variation (June 
and September). The source of the information is the ATLAS OF SURFACE 
CURRENTS /5/. 

Such random behavior lends itself to a probablistic formulation 
rather than a deterministic approach. However the metamorphosis of the 
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frequency plots precludes any distribution other than UNIFORM. This 
assumption is the safest one to make when generalizing for large areas 
and over the whole year. 

A bar chart of current persistency derived from Hydrographic records 
and compiled by Mrs. Martha Olson of USNRDL, San Francisco, Calif. is 
enclosed as figure 28. This chart shows the frequency in large ocean 
areas that the most frequent set in smaller 5° X 5° quadrangles occurs 
for three categories of persistency. For example: In 55% of the 5° X 5° 
sub areas of the Northwestern Pacific the most frequent direction of the 
set occurs in less than 25% of the observations. In 45% of the sub areas 
the most frequent direction of the set occurred in 25% to 50% of the ob- 
servations. Table 7 lists the data from which the bar chart was derived. 

2. Drift. 

Figure 29 was compiled by Mrs. Olson from 1961 Hydrographic 
Office information not yet published. It depicts the frequency of drift 
for various mean speeds of the "most frequent current." It shows that 
in over 80% of the areas the drift of the most frequent current is be- 
tween 0.4 and 0.8 knots. This graph refers only to the North Atlantic. 

Table & that follows lists the mean speed and deviation for various 
sectors at various times of the year. In general the deviation is no 
ans than 0.2 knots, and the mean is less than the mean of the speed in 
the most frequent direction (because when the set is in other than the 
most frequent direction there is a tendency for the movement to be less 
defined. 

C. Conclusions. im 

1. SET is an unpredictable variable at any given time. The safest 

and broadest assumption that can be made is that the direction of current 
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at any given time and at any given place is UNIFORMLY distributed. 

2. The strip considered in figures 20 through 27 was sampled for 
the month of January to determine the percent of area under the curve 
within the semicircle centered at the prevailing set. Two thirds of the 
total area is included under this half of the curve. 

3 No matter what the direction of the set is, the drift is 
reasonably predieteble. The variance of the drift is small and in view 
of the time late involved in the problem, the solution is not sensitive 
to such small variances. 
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TABLE 7 
A MEASURE OF CURRENT PERSISTENCY IN FIVE SELECTED AREAS 
PART 1 


Number of Number of Number of Quadrangle-Months 
Quadrangle- Observations in Which Most Frequent Di- 

REGION Months rection is Observed the 
Indicated % of Observations 


a 25% 252-50% 50% 
A. Northwestern 


Pacific Ocean 80 9,512 45 35 0 
B. Northeastern 

Pacific Ocean 150 12,595 62 78 10 
C. North Atlantic 

Ocean | 253 184, 760 90 . 149 14 
D. East Indian | 

Ocean 59 1,092 27 31 1 
E. Northwestern 

Indian Ocean 96 10,236 28 64 te 

TOTAL 638 218,195 252 357 29 

PART 2 
% of Quadrangle-Mohths in Which Most Frequent Direc- 
REGION tion is Observed the Indicated % of Observations 
2 2 

A. Northwestern 

Pacific Ocean 56% LL% OF 
B. Northeastern 

Pacific Ocean Lid 52% 7% 
GC. North Atlantic 

Ocean 35% 59% 6% 
D. East Indian 

Ocean 66% 29% 5% 
E. Northwestern 

Indian Ocean 29% 67% L% 

TOTAL 39 56 2 


*Additional information on the following page. 
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TABIE 7 — Cont'd. 

Quadrangle-Month: Grouping of data from all observations over the 
years taken during the indicated calendar month and in the same quad- 
rangle of approximately 5°X5°. 

Source of Data: Measurement of lengths and widths of spokes of 
current roses in HO Publications, Atlas of Surface Currents. 


Courtesy of Mrs. Martha Olson, 
USNRDL, San Francisco, Calif. 
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APPENDIX V 
SEARCH PLAN FOR SMALL FIX DEVIATION 
Two other errors in location are considered besides the fix error. 
a. Navigational errors of an aircraft in finding the exact 
location to which dispatched are characterized by a normal distribution. 
In this model it is assumed that the error due to navigation is normally 


distributed with a standard deviation of five miles. The variances of 


the fix error and navigation error are additive with a resulting standard 


deviation of 5.6 miles (}(2.5)*+(5)*). 

be Current flow is assumed to be equally likely in direction 
at a given location and time. The drift is dependable enough as pub- 
lished for the area of interest. 

The first pass of the area should conform to the pattern shown in 
figure 31. The search is designed to investigate the peak of the dis— 
tribution at the time of arrival and follow up with loops outside al 
then inside. This will maximize the probability of detection per time 
of search effort. Koopman /2/ discusses this plan in detail. 


u = drift 


v = aircraft velocity 
T = time elapsed since initial detection 


The time to go from A to B in Figure 30 is r,+tro 
r 
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| 
‘ 


Fig. 30 Pattern for retiring search. 





The time for the object of the search to move from circle r, to 


circle ro = ro-r4 
u 


To make this equal to search vehicle's time to go from A to B, 


Pro = mr,, where m= yiue Similarly rg = mro; ry = mr33 etc.... 
Vou 


To prevent overlap and expand or contract about the loop rs # mr,. 
To make this leg be S miles outside the first leg determine rs such that 


u 


retry, = re-—r,-S, from which rs = mr,+a where a “S 


The distance ry should equal ul. However since in practice the ex=- 


panding square is used to approximate a spiral, a mean r, is selected so 
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that ry = 0.%it. To further allow for the fact that any change in course 
of the pool, at the same speed will reduce the outward component, rj, is 
selected as 0.8ut. 

= mr, 

= my try 
nlp 
nL, +a (outside) 


L, 
Ly 
ly 
ly 
Le = ml, 
ly 
ly 
La 
Ig 


= nigta 
= my 
=mly-2a (inside) 
== 
lao = mig-2a 
at) mag | 
Lo = mLj4+3a (outside ) 
etc... 


A sample search is worked out for the following situation: 


Height of aircraft = 500 feet 
Spacing (S) = 9.61 miles 
Drift (u) = 0.5 knots 


Speed of aircraft (v) = 180 knots 


Time late (T) = 24 hours 





ry = .8(.5)(24) = 9.6 miles 
a = ys = 180(9,61) = 9.65 
v-u 179.5 
% = ta = 180.5 = 1.0056 
v-u 179.5 | 

Iy = 9.65 Total mileage for this pattern 460 mi. 
In = 19.3 Time for pattern ~- 2 hra 35 min. 
In = 19.4 If this plane only has 3 hrs on station 
Ly, = 29.2 available it should spend the last 25 
Is = 29.4 minutes sweeping the lower left corner. 
Le = 39.3 Another aircraft will have to make the 
Ly = 39.5 parallel sweep.secom pass which re- 
Lg = 20.4 quires 1 hr 40 min (300 miles). 
Ig = 20.5 
Lio = 1.3 
Log = 14331 
Lyp = 30.3 
L,3 = 30.5 
In, = 59.7 
ms = 9° 
Lyg¢ (end) = 60 


The plan leaves a holiday in the lower left corner of the pattern. 
In order to take advantage of the best information the upper right corner 
should be oriented in the direction of the prevailing current so that the 
starting leg is 45° to the left of the most frequent current. 

Note that the center may be adequately covered before outside areas 


are searched —- depending on the balance of the time late and standard 
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deviation of location. In case this happens the pattern should be con- 
tinued with outside expansion only. In the example cited above the cen- 
ter search is just a narrow turn of 270° to right. 

After this search is completed proceed with a parallel sweep as ex- 


plained in paragraph 11b of the main body of this study. 
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APPENDIX VI 
SEARCH PLAN FOR LARGE FIX DEVIATION 
For the case where the fix is not any more certain than within a 
standard deviation of 25 miles or more, the errors in navigation and 
current movement can be ignored. A set of expanding squares —- squares 
of uniform coverage - are superimposed in the square pattern. Koopman 
/2/ covers the subject in detail. 
The altitude of the aircraft should be as low as operationally 
acceptable and for an example 500 feet will be used. 
Figure 32 is a drawing of the pattern for the following problem: 
Aircraft speed = 180 kts. 
Endurance on station 2 540 miles 
Spacing = 11.0 miles 
Average search width (W=E) = 3.44 miles 
Coverage factor @ = WL = 4558 square miles 
Maximm search radius = a = (43.7 miles = = #) 


Half length of Kth square (S,) = 4% W G2 ke | 


Number of passes n = a“S_ = 2 
. 25°W 


S, = 19.5 
Sy = 33.8 

Length of a side of K*th square = 25, 
ly = 39 
Ly = 67.6 
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Number of circuits within kth major square (N. ) = Ik 
29 


Ny as 

Ny = 3 

The second uniform square is tilted 45° to the first. This pattern 
is plotted in figure 32. 


The average probability of detection iss: 


Pi 1(1 nll exp (ai = 434 


If the search aircraft has more endurance than considered above (3 
hours at 180 kts - assuming the search area up to 1000 miles from base) 
the plan should be appropriately recomputed. 

For instance if search area were only 500 miles from base 

Station Endurance 1500 
Coverage Factor § = WL = 12660 
Maximum search area radius a = 55.4 


Number of passes n = 3 


54 = We In = 39 N, =2 
So = 33.8 In = 6726 Np = 3 
53 = 43.6 L, = 87.2 Ny =4 


The average probability in this case would be: 


P = 1~(14t ep (ail |= 0.67 
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